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Abstract. I review some results of estimation of the power spectrum of 
density fluctuations from galaxy redshift surveys and discuss advances that 
may be possible with the Sloan Digital Sky Survey. I then examine the re- 
alities of power spectrum estimation in the presence of Galactic extinction, 
photometric errors, galaxy evolution, clustering evolution, and uncertainty 
about the background cosmology. 



1. INTRODUCTION 



The advent of deep redshift surveys of 10 4 — 10 6 galaxies, such as the 
Las Campanas Redshift Survey (LCRS; Shectman et al. 1996) and the AAT 
2df survey (Colless 1998), and multiband photometric and spectroscopic 
surveys such as the Sloan Digital Sky Survey (SDSS; Gunn & Weinberg 
1995; Szalay 1998), marks the beginning of a new era of investigations in 
large-scale structure. Rather than treat galaxies as indistinguishable tracers 
of mass in a static distribution, we will study the detailed dependence of 
clustering on galaxy type and cosmic epoch. In fact, as I will discuss, we 
must study the species and redshift dependence of clustering to precisely 
differentiate among cosmological models. In somewhat shallower wide-angle 
surveys (Huchra et al. 1983; Giovanelli & Haynes 1984; Geller & Huchra 
1989; da Costa et al. 1994; Fisher et al. 1995) it has been possible to begin 
study of these effects; at redshift z = 0.1 and larger it becomes critical to do 
so. In addition to precisely characterizing galaxy clustering at the present 
epoch, we will use the anisotropy of clustering, as well as the evolution of 
this clustering, to constrain cosmological parameters. 
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A key issue that I examine is the limiting precision with which we can 
hope to measure galaxy clustering at the present epoch. Although much has 
been made of the power of large galaxy surveys to measure cosmological 
parameters, most of this discussion has focussed on finding an analysis 
method that gives "optimal" uncertainties on model parameters, where 
the uncertainties are due to cosmic variance and shot noise (Vogeley & 
Szalay 1996; Tegmark et al. 1998). Much less has been said about the 
uncertainties caused by random and systematic errors in photometry or 
Galactic extinction corrections, and systematic uncertainties that arise from 
issues of cosmology, and evolution of galaxies and their clustering. In this 
paper I discuss the impact of these effects on analysis of galaxy redshift 
samples from the SDSS. 

2. MEASUREMENTS OF THE POWER SPECTRUM FROM 
GALAXY SURVEYS 

Rather than give a review of the myriad measures of galaxy clustering 
that we hope to apply to these new surveys, here I focus on one statistic 
and examine how possible random and systematic uncertainties will affect 
the accuracy of this measurement. Not surprisingly, I consider the power 
spectrum of galaxy density fluctuations, because lowest-order measures are 
where we first gain precision, and because uncertainties in higher-order 
measures often scale from the uncertainties in the variance. 

2.1. THE DATA 

Figure 1 plots the power spectrum P(k) in units of h~ 3 Mpc 3 (a power 
spectral density) as estimated from a number of radio (Peacock & Nichol- 
son 1991), infrared (Fisher et al. 1993; Tadros & Efstathiou 1995), and 
optically-selected (da Costa et al. 1994; Lin et al. 1995) galaxy redshift sur- 
vey samples, and the real-space power spectrum inferred from angular clus- 
tering of a photometric sample (Baugh & Efstathiou 1993). To avoid hope- 
lessly confusing the plot, no error bars are shown (unfortunately, the various 
authors differ somewhat in their choice of binning and methods of comput- 
ing these uncertainties, so comparison of these errors would not be as in- 
formative as one might hope). The shapes of most of these spectra roughly 
agree over the decade 10 — lOO/i -1 Mpc. However, there is large variation 
of their amplitudes. We usually attribute this variation to differences in the 
"bias parameter," which affects the amplitudes as P\(k)/P2{k) = (&1/62) 2 
if the bias is scale-independent. This strong dependence of clustering am- 
plitude on method of galaxy selection implies that detailed knowledge of 
the selection criteria and construction of sub-samples for analysis that are 
homogeneous in this selection are prerequisites for precision measurement 
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Figure 1. Estimates of the redshift-space power spectrum from a variety of redshift sur- 
veys, and an estimate of the real-space power spectrum inferred from angular correlations 
(APM). 



of the power spectrum. 

To clarify comparison of the shapes of the power spectra, Figure 2 
plots these same data, with the amplitudes of the curves shifted so that 
all match P{k) of the LCRS sample near wavenumber k = 0.1. No shift 
was required to obtain agreement between the LCRS and the CfA2+SSRS2- 
100 volume-limited sample, which is not surprising since they both include 
optically-selected galaxies that are roughly M* and brighter. In general, one 
sees excellent agreement between optical redshift samples for wavenumbers 
k > 0.1. There is some disagreement between infrared and optical samples 
over this same range of wavenumber, perhaps because the selected galax- 
ies sample different physical environments. One should be careful not to 
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Figure 2. Power spectrum estimates from Figure 1 arbitrarily scaled to match at k = 0.1 
to allow comparison of power spectrum shapes. CfA2+SSRS2-101 and LCRS were not 
adjusted. 



over-interpret this comparison of spectral shapes. Where large shifts are 
necessary to match the amplitudes, this shifting procedure is questionable, 
because the effects of redshift distortions (both the linear boost described by 
Kaiser 1987 and the washing out of small-scale power due to redshift fingers 
of clusters) and non-linear evolution both depend on the ratio of galaxy- 
clustering to mass-clustering amplitude. On large scales, at wavenumber 
k < 0.1, we find significant departures in the shapes of the spectra, with a 
scatter of roughly a factor of three among the samples. This scatter is con- 
sistent with the uncertainty due to cosmic variance (small sample volume 
relative to the wavelength scales of interest). 
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Figure 3. Non-linear redshift-space power spectra for a variety of COBE-normalized 
CDM models. 



2.2. MODEL FITS 

The shape of these spectra may be well-fitted by a linear Cold Dark Mat- 
ter (CDM) power spectrum (Bardeen et al. 1986) with shape parameter 
T = 0.25 (where V = Qh in the simplest models, with h = H /W0). Here 
we use the linear CDM spectrum merely as a fitting function. It is inter- 
esting to note that, for values of T = flh ~ 0.25, this linear curve does a 
reasonable job in fitting the redshift-space power spectrum of a fully non- 
linear evolved CDM model because the effects of non-linear evolution and 
redshift distortion nearly cancel out. 

Does this fit to the spectral shape imply that our work is done? Hardly. 
This rough shape might be matched by tweaking the parameters of several 
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Figure 4- Real-space power spectra for the same models plotted in Figure 3. 



competing models (e.g., by adding some cosmic density in neutrinos and/or 
by varying the slope of the primordial power spectrum). We require higher 
resolution to detect features in the power spectrum that might be caused 
by acoustic oscillations near recombination (Eisenstein & Hu 1998) or other 
physical effects. Higher resolution in the Fourier domain requires a larger 
survey volume (see next section), which is also necessary to probe large 
wavelength scales. Measurements of the power spectrum to accuracies of a 
few percent on scales A > 100/i _1 Mpc would allow detailed comparison of 
clustering of galaxies at the present epoch with clustering of the mass at 
redshift z = 10 3 , as revealed by the anisotropy of the Cosmic Microwave 
Background (CMB). 

Figure 3 illustrates how CDM models with different parameters yield 



HIGH-PRECISION LARGE-SCALE STRUCTURE 



7 



redshift-space power spectra that are similar. To compute these spectra, I 
use the semi-analytic formalism described by Peacock (1997) to approxi- 
mate the effects of non-linear gravitational evolution and redshift distor- 
tions on the linear CDM spectrum. All of these models have a mass spec- 
trum that is normalized to match the COBE DMR measurement of the 
CMB anisotropy (Wright et al. 1996). Even without exploring admixtures 
of CDM and neutrinos (e.g., Primack et al. 1995), we see that several choices 
of parameters might fit the equally well (or equally poorly, depending on 
one's taste). Not examined here is how the small-scale power spectrum de- 
pends on the details of galaxy formation relative to the mass distribution. 

2.3. REAL VS. REDSHIFT SPACE 

One of the problems of working in redshift space is that a large amount of 
small-scale power creates a large small-scale velocity dispersion. In redshift 
space, the resulting "fingers of God" wash out much of this same small- 
scale power. Similarly, for fixed normalization and Hubble constant, CDM 
models with large £1 experience a significant boost of their redshift-space 
power spectrum amplitudes (1.87 if 0, = 1 and b = 1), which makes their 
large-scale spectra look more like those of low Q models. 

Redshift-space power spectrum measurements with accuracy of a few 
percent could differentiate between the model spectra in Figure 3 but, 
clearly, it would be easier to break the model degeneracy that we see in 
redshift-space power spectra if we could access the underlying real-space 
(i.e., configuration space) power spectrum of the galaxy distribution. Fig- 
ure 4 shows the power spectra of the same models from Figure 3, now 
plotted in real space. 

The lesson here is not that we should use the angular correlation func- 
tion, and give up measuring redshifts, but rather that we should take advan- 
tage of the anisotropy of clustering in redshift space to allow measurement 
of both the power spectrum and 0. Andrew Hamilton's contribution to this 
volume examines this issue in great detail. 

3. MINIMAL POWER SPECTRUM UNCERTAINTIES FOR 
THE SDSS 

Power spectrum measurements with accuracy of a few percent on scales 
from one to a few hundred Mpc would allow us to make great strides in 
constraining cosmological models. Is this feasible? Here I begin to examine 
the accuracy with which we dare hope to measure P(k) in the foreseeable 
future. 
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Figure 5. Predicted uncertainties in the power spectrum estimated from a vol- 
ume-limited (Rmax = 500ft _1 Mpc) sample of SDSS North and for the Bright Red Galaxy 
sample (upper set of error bars). These errors assume that the true power spectrum is 
that of an Slh = 0.25 CDM model and that the BRGs are more clustered than normal 
galaxies. Plotted for comparison to the SDSS North errors are CDM power spectra (nor- 
malized to cr 8 = 1) for a range of Q,h from 0.2 (uppermost curve) to 0.5 (lowest curve) 
and indications of the range of comoving scales probed by the COBE and MAP CMB 
anisotropy experiments. 



3.1. APPROXIMATE UNCERTAINTY PREDICTIONS 



For fixed survey strategy, the survey volume and the number density of 
galaxies in the redshift sample set a lower bound on the uncertainty in 
the estimated power spectrum. In the limit of a perfectly spherical volume- 
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limited sample, the uncertainty in the estimated power per mode is roughly 

S(k) 




1 + 



P{k) 



(1) 



(Feldman, Kaiser, & Peacock 1994). P(k) and P(k) are the true and es- 
timated power spectra, respectively. S(k) = 1/n is the shot noise power 
for mean galaxy density n. V c = (2tt) 3 /Vs is the coherence volume in the 
Fourier domain for a survey with volume Vs- It is assumed that we average 
the power estimates over a shell in Fourier space with volume 47rk 2 Ak. 
That is, we average the power over all angles, and over bins with width 
Ak > 2tt/R, where R is the survey depth. For a larger volume, the coher- 
ence volume in Fourier space decreases, thus we obtain a larger number of 
independent probes of the power spectrum. This convenient approximation 
breaks down when the survey volume is highly anisotropic; in this case 
the coherence volume in Fourier space (the "window function") becomes 
anisotropic and the range of true wavenumber that is probed depends on 
the direction of k. Figure 8 of Alex Szalay's contribution to this volume 
shows window function shapes and relative sizes for a variety of redshift 
surveys, including the SDSS. 

A volume-limited sub-sample of the SDSS redshift survey in the North 
Galactic Cap region (this covers solid angle of n steradians) that has co- 
moving coordinate depth of R = 500/i -1 Mpc will include all galaxies with 
absolute magnitude roughly 0.4 m brighter than M* , with mean number 
density n ~ 2 x lCT 3 /i 3 Mpc~ 3 . Hereafter I will use this as a fiducial sample 
for computing uncertainties in the estimated power spectrum. This sample 
is conservative, in the sense that the SDSS will allow selection of much 
deeper samples, which include more volume, and because a magnitude- 
limited sample to the same depth would have larger average number den- 
sity, hence smaller shot noise. 

Figure 5 plots la uncertainties for power estimates binned by Ak = 
27r/(500/i~ 1 Mpc) and assuming that the true power spectrum is described 
by a T = 0.25 linear CDM spectrum. Also shown are a family of CDM 
spectra with different shape parameter T. The uncertainty in each bin of 
width Ak at a wavelength scale of A = lOO/i^Mpc will be a(P)/P « 0.06. 
Such a measurement would easily differentiate between various CDM mod- 
els, as the family of CDM curves indicates. Comparison with results from 
the LCRS further frustrates the power of such a large survey: if a "bump" 
at A ~ 128/t _1 Mpc were found in the SDSS sample with SP / 'P sm ooth = 0.76 
in a single bin (as seen by Landy et al. 1996), this would be a 12<r event. 
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3.2. SAMPLES WITH LARGER VOLUME? 

Could we do better than the SDSS with, for example, an all-sky survey to 
similar depth? This suggestion is not out of the question, if we combine all 
of the redshift surveys that will exist in a few years' time. The SDSS NGP 
survey will cover one quarter of the sky; a similar full-sky survey would have 
power spectrum uncertainties that are at best a factor of two smaller. More 
realistic at optical wavelengths would be a survey over Galactic latitude 
\b\ > 30, which covers half the sky, yielding errors that are at best 40% 
smaller than the SDSS. In fact, the SDSS itself will cover an effective area 
larger than ir steradians, by including three slices in the southern Galactic 
hemisphere, each 3°x ~ 100° in area. If the SDSS is combined with the 
AAT 2df redshift survey and the LCRS, then most of the high-latitude sky 
will be sampled. Galaxy selection in the infrared would allow a true full- 
sky survey, but no IR photometric survey exists to depth comparable to 
the SDSS. The IRAS PSC-Z survey (Saunders et al. 1994) comes closest 
to this goal, but has much larger shot noise (smaller galaxy density) and 
apparently more severe sample evolution. 

One-tenth of the million galaxies targeted for spectroscopy by the SDSS 
will be a sample of Bright Red Galaxies (BRGs) that are volume-limited 
to redshift z = 0.45. The volume limit will be enforced by selecting the 
galaxies by their absolute magnitudes, estimated using the photometric 
redshift technique (Connolly et al. 1995). This sample will cover a volume 
eight times larger than the fiducial SDSS-500 sample described above, albeit 
with larger shot noise. The upper curve and error bars in Figure 5 show the 
expected uncertainties from cosmic variance and shot noise in the SDSS 
BRG sample, using equation 1 and assuming that the BRGs are biased by 
a factor of 2 (factor of 4 in the power spectrum) relative to the normal 
galaxy sample. The brightest and reddest galaxies appear to evolve least of 
all, so this sample will provide our best chance to measure the evolution of 
clustering (see section 5.3 below). 

4. MEASUREMENT ERRORS AND EXTERNAL SYSTEM- 
ATIC S 

In the previous section I describe "ideal" uncertainties that are lower bounds 
on the total error budget. I now examine whether these expectations are re- 
alistic. What other effects limit our knowledge of clustering at the present 
epoch? I use the SDSS as a test case, but the following analysis is also 
relevant for other surveys to similar depth. 

There are several possible effects about which we do not need to worry 
for the SDSS. Because the galaxies in the spectroscopic sample will be 
five magnitudes brighter than the point source detection limit, star-galaxy 
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separation will not be a problem. Galaxies with very low central surface 
brightness will be excluded from the redshift sample, because we would not 
be able to get redshifts of these objects, but we will understand this se- 
lection bias quite well. For example, our magnitude and surface brightness 
cuts will be "fuzzy," such that we observe some fraction of objects beyond 
the nominal cuts. The redshifts that we do measure will have 20 — 30 km s _1 
accuracy, so this project is overkill for simply measuring redshifts and we 
expect very few failures. The use of a fiber-fed spectrograph with plugplates 
to hold the fibers in the focal plane imposes a 55 arcsecond minimum sepa- 
ration for pairs of objects, except where a target lies in the overlap between 
two or more plugplate fields. However, during the Fall observing season, 
the larger ratio of spectroscopic time to photometric time (because we will 
image a smaller area), will allow a larger covering factor for spectroscopy, 
thus we will obtain spectra for at least all pairs of objects and potentially 
could observe all n-tuples up to n = 7. From these southern data, we will 
study the effects of excluding close pairs from the northern spectroscopy 
(e.g., we can study the very small-scale velocity dispersion of galaxies from 
the southern spectroscopy). 

4.1. EFFECTS OF EXTINCTION AND PHOTOMETRIC ERRORS 

Despite our best efforts at correcting for Galactic extinction and calibrating 
the photometry, errors in these will add to the minimal uncertainties illus- 
trated in Figure 5. In collaboration with Andy Connolly at Johns Hopkins, 
I have investigated the impact of such errors on power spectrum estimation 
from samples of the SDSS (Vogeley & Connolly 1998). The remainder of 
this section is an overview of some of our results. 

Errors in either photometry or the extinction corrections affect the ap- 
parent number density of galaxies in similar fashion. Here we consider only 
the effects of extinction or photometric error on total magnitude selection. 
Although the SDSS spectroscopic samples will also be chosen on the ba- 
sis of a central surface brightness cut, simulations of the SDSS indicate 
that surface brightness dimming from extinction will affect an insignificant 
number of galaxies. 

In a volume-limited sample that should include all galaxies brighter 
than absolute magnitude Mn m , a magnitude error of Am over a patch on 
the sky causes the apparent galaxy density in that direction to differ from 
the true galaxy density by 



n obs 



$[< M - Am] 



(2) 



n true 



$[< M] 
(p(M)Am 
$[< M] 
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Figure 6. Effects on the observed galaxy density that would result from 0.1 m of ex- 
tinction in a magnitude-limited (upper panel) or volume-limited (lower panel) redshift 
sample of the SDSS. 



~ exp -AmlU"^" " ±lim > , 
where <j)(M) and <E>(< M) are the differential and integrated luminosity 
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functions of the galaxies. The approximations are valid for a Schechter lu- 
minosity function when Mu m is close to or brighter than M*. This form 
clearly shows how the strength of the effect depends on the slope of the 
luminosity function near the magnitude limit; deeper sub-samples will be 
more severely affected by extinction or photometric uncertainty than shal- 
low samples because the count slope is steeper. For a magnitude-limited 
sample with limit mn m , the absolute magnitude limit varies with distance 
Mu m (r) and the apparent density differs from the true galaxy density by 

n obs = $[< M(r) - Am] 

n true *[< M(r)] • ( ' 

Figure 6 shows the effect of a Am = 0.1 m error on the apparent density 
of volume-limited and magnitude-limited samples of the SDSS. This mod- 
ulation of the apparent density grows with distance in a magnitude-limited 
sample and is, on average, less severe than in an volume-limited sample that 
is cut off at the same depth. In the following analysis I consider volume- 
limited samples because (1) the effects are more severe for these, thus the 
analysis is conservative in the sense of giving worst-case errors and (2) we 
won't believe the results of magnitude-limited power spectrum analyses un- 
til we understand how the clustering amplitude differs among galaxy species 
(see section 2.1 above). 

4.2. THE POWER SPECTRUM OF GALACTIC EXTINCTION 

Because Galactic extinction modulates the apparent density of galaxies, 
it could cause erroneous clustering power to appear in our data. To lead- 
ing order, the net effect is to add to the apparent clustering, P b s (k) ~ 
Ptrue(k) + Pextinction(k) ■ If we make no correction at all for Galactic ex- 
tinction, we expect to observe extra clustering power in the SDSS-500 sam- 
ple, as shown in Figure 7. Here we predict the effect of extinction over the 
SDSS North region using a cosecant law, the Burstein & Heiles (1982) map, 
and the Stark et al. (1992) HI map. These maps predict similar clustering 
power; it is the large-scale cosecant-like variation in extinction that causes 
the trouble. Restricting the analysis to somewhat higher Galactic latitude 
(|6| > 40°) only slightly ameliorates the problem. Failure to notice that we 
live in a Galaxy that is laced with dust would cause one to infer that the 
power spectrum of the universe suddenly rises in power-law fashion beyond 
200 - 300k- 1 Mpc. 

Fortunately, we are not oblivious to the effects of dust and we plan to 
make some correction for this effect. The SDSS will select spectroscopic 
targets after applying an a priori correction for extinction to the apparent 
magnitudes, thus constructing a sample that is uniform in an extragalac- 
tic sense. As of this date we plan to use the extinction maps that have 
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Figure 7. False observed clustering power that would be caused by Galactic extinction, 
for different assumed extinction maps. For comparison is an Qh = 0.25 CDM power 
spectrum and the expected uncertainties for a volumedimited sample of SDSS North. If 
we make no correction for extinction, then the false clustering power exceeds the true 
power for A > 300fe _1 Mpc. Cutting back the edges of the survey to b > 40° has only a 
small effect. If we correct the B&H map using a crudely-calibrated map derived from the 
Stark et al. HI maps, the residual power exceeds the true power for scales A > 700/i _1 
Mpc. 

been constructed from the DIRBE and IRAS satellite data by Schlegel, 
Finkbeiner, & Davis (1998). If we or others later construct a better ex- 
tinction map, we will take the residuals into account when we analyze the 
data. The erroneous clustering power that we really need to worry about 
is caused by unknown residuals between our best extinction map and the 
true Galactic extinction. 

The lower dot-dashed line in Figure 7 shows what happens if the true 
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Figure 8. False clustering power due to fluctuations of the photometric zeropoint (or 
errors in the extinction map) on 3° and 6° scales, or between stripes of the SDSS. Also 
shown is the effect of a 0.1 m gradient across the sky. 

extinction is described by the Burstein-Heiles map and we use a crudely- 
calibrated (here we applied our own ad hoc dust-to-gas relation) version of 
the Stark et al. HI map to estimate and correct for the extinction. Even with 
this roughly calibrated map, the extinction power is below the true power 
on scales smaller than 700/i -1 Mpc. A well-calibrated map would allow us 
to accurately probe the entire range of accessible scales. 

4.3. THE POWER SPECTRUM OF PHOTOMETRIC ERRORS 

Similar to the effects of extinction, errors in photometry that are correlated 
over a patch of sky will modulate the apparent density of galaxies. Ran- 
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dom (uncorrelated between galaxies) errors will contribute to the notorious 
Malmquist bias; here we consider only correlated photometric errors, such 
as those that might arise from variations in the zero point. 

Figure 8 shows the extra clustering power in the SDSS-500 sample that 
we would observe if the photometric zero-point varies by 0.1 m over patches 
that are 3° in diameter (the field of view of the SDSS camera), or by 0.05 m 
over patches that are 6° in diameter (the size of a Schmidt plate). These 
curves should match on scales much larger than the patch size, but vary 
somewhat because we compute these curves using a Monte Carlo realiza- 
tion of the errors. Even such gross errors in photometry produce erroneous 
clustering that is less than 10% of the true power for A < 10 3 /i _1 Mpc. 

The relatively small amplitude of power caused by even 0.1 m photome- 
try errors over the field of view of the SDSS camera is very good news for 
the extinction problem. This result implies that we can cure the extinction 
problem if we can use the SDSS data themselves to construct an extinction 
map that has random field-field errors that are smaller than 0.1 m (but, of 
course, without any large-scale systematic errors). Variations in faint galaxy 
counts and colors and colors of hot Galactic halo subdwarfs offer several 
means for constructing such a map, which would provide an independent 
method for computing extinction corrections. 

The lesson of the extinction errors is that large-scale gradients, rather 
than small-scale random photometry errors, are the real worry. Figure 8 also 
shows the extra clustering power that arises if there is a zero-point gradient 
across the sky (pole to pole) of 0.1 m or if the zero-point randomly varies 
by 0.05 m between the "stripes" of the SDSS. In other words, the latter 
illustrates what happens if the photometry is consistent within each ~ 
3° x 120° scan region, but the calibration zero-point varies between disjoint 
regions. Even such a severe error (which we would no doubt notice by 
comparing the overlap regions of the scans) would produce extra power 
that is below 10% of the true power on scales < 10 3 ft, _1 Mpc. To be certain, 
we would not be happy to have a 10% effect in our measurements. The 
point is that this is clearly an upper bound on the magnitude of such an 
effect. 

It is important to note that, just as examination of the anisotropy of 
clustering will reveal the effects of redshift distortion on the power spec- 
trum, comparison of clustering in the angular and radial directions will 
provide a test of false clustering due to extinction or photometric error. 
Cross-correlation of the angular clustering of samples that are selected to 
lie at different distance would clearly reveal the signature of correlated 
magnitude errors. 
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5. EVOLUTION: THE UNIVERSE, GALAXIES, AND CLUS- 
TERING 

What is the "present epoch?" For statistical analyses of clustering, a useful 
definition is the redshift range beyond which ignorance of, or inability to 
correct for, the effects of evolution of galaxies and their clustering becomes 
the dominant source of uncertainty. In this section I examine how the ap- 
parent power spectrum is affected by (1) the mapping between redshift and 
comoving coordinate distance, (2) galaxy evolution, and (3) clustering evo- 
lution. Careful analysis will allow us to learn about evolution; the question 
here is, when do these effects begin to dominate the apparent clustering? 



5.1. KNOW THY COSMOLOGY 

The relationship between redshift and comoving coordinate distance de- 
pends on the deceleration parameter qo, 

r ( z ) = M -27] T 3 W z + (! - 9o) (l " a/1 + 2(Zo*)] • (4) 
-"o^oli + z ) 1 v ' J 

Out to redshifts of a few thousand km s _1 , it matters very little which value 
of qo we use. But at larger distance, this transformation begins to affect 
galaxy clustering in several ways. Using the wrong qo will cause fluctuations 
on a fixed comoving scale to appear as fluctuations on different apparent 
scales as a function of redshift. Over the redshift range probed by the SDSS 
this effect would cause only minor smearing of features in the spectrum; 
varying qo from 0.1 to 0.5 changes the length scale by only 4% over the 
redshift range z = to z = 0.1. 

It also follows that the comoving volume element depends on qo, as 

= 2vr (±\ 3 [go* + (l-go)(l-\/l + 2goz)] 2 , 

If we assume a value of qo that is too large, the volume element dV will 
erroneously decrease with distance, and cause a radial gradient in the ap- 
parent galaxy density. Varying qo from 0.1 to 0.5 raises the comoving galaxy 
density by 8% at redshift z = 0.1. 

Because construction of volume-limited samples relies on using the red- 
shift to compute the absolute magnitudes of the galaxies, an error in qo 
translates into an error in the inferred absolute magnitudes, 

AM = 5 log \r(z, ql rue )/r(z, gg Mwmed )l . (6) 



This error in the absolute magnitudes yields a gradient in apparent galaxy 
density, as described by equation 2. Using too large a value of qo, we would 
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Figure 9. Density gradients caused by (from top to bottom) galaxy luminosity evolution 
with ft — 2.3, density evolution with a — 2, error in the computed comoving volume, and 
error in the computed absolute magnitudes in an absolute-magnitude limited sample of 
the SDSS. The latter two cosmological effects are shown for the case in which the true 
go = 0.1 but we have used qo = 0.5 to relate redshifts to comoving coordinate distances. 



infer that galaxies are fainter than their true luminosity and exclude some 
from our sample, thus lowering the galaxy density at large redshift. At 
redshift z = 0.1, varying go from 0.1 to 0.5 lowers the density by 7%. 

Figure 9 shows the density gradients in the SDSS-500 volume-limited 
sample that would result from assuming qo = 0.5 if the true value is qo = 
0.1. Interestingly, the rise in apparent density that results from the change 
in comoving volume element nearly cancels the absolute-magnitude effect 
(this cancellation is not universally true because the absolute-magnitude 
effect depends on the slope of the integrated luminosity function at the 
absolute magnitude limit). 
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k [2tt/X] 



Figure 10. False clustering power that would be observed in the SDSS-500 vol- 
ume-limited sample by (top to bottom) galaxy density evolution with a = 7, galaxy 
luminosity evolution with (3 = 2.3, density evolution with a = 2, and from gradients 
caused by using the wrong value of go- In every case, the sharp turnover at small k is an 
artifact of our failure to correct for a numerical effect. 



Figure 10 illustrates the extra clustering power that results from these 
effects, considered separately. Both the comoving volume effect and the 
absolute magnitude effect cause negligible power because the gradient is so 
shallow. Even if these gradients did not cancel, they would not be a cause 
for worry in samples that extend to R = 500/i _1 Mpc. 

5.2. GALAXY EVOLUTION 

Also shown in Figures 9 and 10 are the density gradients and extra clus- 
tering power predicted for galaxy luminosity or density evolution. Here we 
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model density evolution as n(z) oc (1 + z) a with a = 2 shown in Figure 9, 
and both a = 7 (upper curve) and a = 2 (middle curve) in Figure 10. If 
a is much larger than 2, this effect could cause trouble in the SDSS sam- 
ples. For example, Saunders et al. (1990) find evidence for a = 7 in the 
IRAS-selected galaxy luminosity function. We model luminosity evolution 
as L(z) oc (1 + z)P with (3 = 2.3, as fit by Lilly et al. (1995) to their 
/-band luminosity functions. Saunders et al. find (3 = 3 for IRAS galaxies. 
Figure 10 shows that this type of luminosity evolution would contribute a 
significant amount of clustering power on the largest scales probed by the 
SDSS-500 sample. 

A related effect is the difficulty of computing the K-corrections for galax- 
ies over a large range of redshift. Out to redshift z = 0.2, failing to perform 
the proper K-corrections in r' (the selection band for SDSS spectroscopy) 
causes an effect on the galaxy samples that is somewhat smaller than the 
luminosity evolution effect. 

The lesson here is that one cannot treat galaxies as indistinguishable 
"points" for the purpose of power spectrum analysis. We already know this 
from, for example, luminosity bias, whereby very bright galaxies are more 
strongly clustered than the average (e.g., Park et al. 1994). As we push 
to larger redshift, we must first characterize the redshift dependence of 
the intrinsic properties of galaxies. Only then can we compute with con- 
fidence the clustering statistics of samples that are either homogeneous in 
their clustering properties or that have known dependence of clustering on 
galaxy species. Computing the requisite multivariate luminosity functions 
and spectral evolution of galaxies requires multiband photometry and ho- 
mogeneous samples with good resolution spectra. The SDSS will provide 
both of these desiderata. 

5.3. CLUSTERING EVOLUTION 

Clustering of galaxies fails to be ergodic as we probe deeper into the universe 
for yet another reason: the pattern of clustering itself evolves with time. In 
linear theory, the growing-mode pertubations grow with the scale factor a 
of the universe such that P(k,a) oc D 2 (a), where D(a) is a function of the 
cosmic matter density and cosmological constant f^A- Relative to the 
present epoch, when D(z) = 1, the growth factor when the universe had 
relative scale size a = (1 + z)~ l is (Carroll, Press, & Turner 1992) 



2 a \a J 



nl/2 



D(a) = 



1) 



(7) 




3/2 



x ss 



a 



da'. 
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Figure 11. Evolution of the clustering amplitude with redshift for (top to bottom) 

{a M = i, fi A = o}, {o M = o.4, Ov = o.6}, {n M = 0.4, n A = 0}, {a M = 0.2, n A = 0.8}, 

and {£Im = 0.2, 17a = 0}. Plotted here is the difference in clustering amplitude relative 
to z = 0. 

For fijvf = 1 and S7a = 0, this relation is simply D(z) = (1 + z)^ 1 . 

Figure 1 1 plots the apparent decrement of clustering amplitude relative 
to z = 0, D 2 (z = 1) — D 2 {z) = 1 — D 2 (z), for a few choices of spatially 
flat and open cosmologies. By redshift z = 0.2, the decrement is 30% for 
Um = 1 and 12% for an VLm = 0.2 open universe. The difference in the 
decrement for these models is then larger than the statistical errors in the 
power spectrum on most scales for a SDSS sample to this depth. 

In other words, the cosmology-dependence of the growth factor becomes 
the dominant source of uncertainty in estimation of the power spectrum 
and we reach a fundamental limit on our ability to measure the power 
spectrum in an model-independent fashion. This is not a problem when 
comparing with models that specify the cosmology, because one simply 
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predicts the power spectrum that would be observed over the specified 
volume for each model. But it is not possible to compare local measures 
of clustering with observations at larger redshift without specifying the 
cosmology. For example, apparent variation in the clustering amplitude of 
a galaxy species that are due to galaxy evolution might be degenerate with 
variation of the cosmology. Thus ends the "present epoch" for the purpose 
of statistical large-scale structure. 

On the other hand, given the large divergence of clustering amplitude 
for models that might otherwise match at z = 0, the strong cosmology- 
dependence of the growth factor seems like a powerful way to constrain 
£Im and ^a- Given all of the caveats above regarding galaxy evolution, we 
must find a sample of objects that (1) has negligible (or, at least, very well 
understood) intrinsic evolution out to a redshift of a few tenths (2) are 
sufficiently numerous that shot noise does not dominate over the cosmic 
variance, and for which we (3) can easily obtain redshifts (either spec- 
troscopic or photometric). One must use a single volume-limited sample, 
because the decrease of clustering amplitude with redshift might be masked 
by the effects of luminosity bias (brighter galaxies cluster more strongly) in 
a magnitude-limited sample. A good candidate for such an analysis will be 
the BRG sample of the SDSS, which will extend to z = 0.45. The bright- 
est, reddest galaxies evolve the least and have strong spectral features that 
ease redshift determination for these objects. Dividing this redshift sample 
into broad bins of redshift and averaging the power estimates over broad 
bands should allow determination of the clustering amplitude to within a 
few percent at each redshift. 

6. CONCLUSIONS 

The large sample volume and dense sampling of the galaxy distribution in 
the SDSS redshift survey will set lower bounds to the uncertainties that 
could, in principle, allow detection of small deviations from smooth spectra 
and easily differentiate between similar cosmological models. Clever use of 
the anisotropy of clustering would allow simultaneous estimation of both 
Qm and the real-space power spectrum. Further, we could use the growth 
rate of clustering to constrain fi^ and f^A in a different way. However, 
all of these dreams hinge on our ability to cleanly separate true galaxy 
clustering from contamination by extinction and photometric error, and 
to differentiate between clustering evolution and evolution of the galaxies 
themselves. 

Thus, the future of large-scale structure lies in studying evolution, both 
of galaxies and of the clustering pattern itself. This evolution is strongly 
species-dependent and will require estimation of multivariate luminosity 



HIGH-PRECISION LARGE-SCALE STRUCTURE 



23 



functions to interpret. In a few years' time, the SDSS, together with the 
2MASS infrared survey and FIRST radio surveys, which will also cover 
(really, have already begun to cover) the same region of sky, will provide 
a database of optical, near-IR, and radio information about millions of 
objects and will be ideal for just this sort of investigation. 

The author acknowledges the collaborative effort of many of the par- 
ticipants in the SDSS, in particular Andy Connolly and Alex Szalay, that 
contributed to the results presented in this volume. Support for this work 
was provided by NASA through grant HF-01078.01-94A from the Space 
Telescope Science Institute, which is operated by AURA, Inc. under NASA 
contract NAS5-26555, 
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